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Abstract 

There exists a distinctive class of physically significant evolution PDEs in one 
spatial dimension which can be treated analytically. A prototypical example of 
this class (which is referred to as integrable), is the Korteweg-de Vries equation. 
Integrable PDEs on the line can be analyzed by the so-called Inverse Scattering 
Transform (1ST) method. A particularly powerful aspect of the 1ST is its ability 
to predict the large t behavior of the solution. Namely, starting with initial 
data u{x,0), 1ST implies that the solution u{x,t) asymptotes to a collection of 
solitons as t — ?• oo, x/t = 0(1); moreover the shapes and speeds of these solitons 
can be computed from u{x, 0) using only linear operations. One of the most 
important developments in this area has been the generalization of the 1ST from 
initial to initial-boundary value (IBV) problems formulated on the half-line. It 
can be shown that again u{x, t) asymptotes into a collection of solitons, where 
now the shapes and the speeds of these solitons depend both on u{x, 0) and on 
the given boundary conditions at a; = 0. A major complication of IBV problems 
is that the computation of the shapes and speeds of the solitons involves the 
solution of a nonlinear Volterra integral equation. However, for a certain class of 
boundary conditions, called linearizable, this complication can be bypassed and 
the relevant computation is as effective as in the case of the problem on the line. 
Here, after reviewing the general theory for KdV, we analyze three different types 
of linearizable boundary conditions. For these cases, the initial conditions are: 
(a) restrictions of one and two soliton solutions at i = 0; (b) profiles of certain 
exponential type; (c) box-shaped profiles. For each of these cases, by computing 
explicitly the shapes and the speeds of the asymptotic solitons, we elucidate the 
influence of the boundary. 
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1 Introduction 



Systems governed by nonlinear differential equations are of fundamental importance 
in all branches of science. Although our understanding of nonlinear equations in 
general is still very limited, the last few decades have seen some important develop- 
ments in the ways we describe, model, and predict certain nonlinear phenomena. One 
of the most exciting developments in this direction has been the implementation of 
the Inverse Scattering Transform (1ST) for integrable nonlinear evolution equations 
in 1 + 1 dimensions and the associated emergence of soliton theory. A particularly 
powerful aspect of the 1ST is its ability to predict the large time asymptotics of the 
solution. Given an integrable evolution equation and initial data uq^x) defined for 
X G M, the 1ST formalism implies that for large t the solution u{x,t) will asymptote 
to a collection of solitons traveling at constant speeds. Moreover, these asymptotic 
solitons can be determined from the initial data: Each soliton corresponds to a zero 
of an analytic spectral function a{k) defined in terms of uo{x) by means of purely 
linear operations. 

One of the most important developments of soliton theory has been the generaliza- 
tion of the 1ST formalism from initial value to initial-boundary value (IBV) problems 
[HIISKTT]. Such problems appear in many applications, where it is often more natural 
to assume that the space variable is defined only on part of the real axis. For exam- 
ple, in the context of water waves it is usually far easier to measure the elevation of 
the water's surface over time at a fixed position (say at x = 0), than to determine 
the surface for all values of x at a specific instant of time. This leads naturally to 
the mathematical formulation of an IBV problem on the half-line {x > 0} with a 
prescribed boundary condition at x = and with vanishing initial conditions. 

For IBV problems, an important issue, both for mathematical and physical con- 
siderations, is the study of the effect of the boundary on the asymptotic behavior of 
the solution. The approach of [H |9] yields an expression for the solution n(x, t) of 
an IBV problem in terms of the solution of a Riemann-Hilbert (RH) problem. All 
the ingredients of this RH problem are defined in terms of the initial and boundary 
values of u{x,t). The asymptotic behavior of u{x,t) can be analyzed in an effective 
way by using this RH problem and by employing a nonlinear version of the steepest 
descent method [7J. Using this approach it can be shown ^31 113 US] that, just 
like in the infinite line setting, the solution on the half-line will asymptote for large t 
to a collection of solitons traveling at constant speeds. It should be emphasized that 
on the half-line the asymptotic solitons are not generated by zeros of the equivalent 
of a{k) (which is specified by the initial conditions), but they are generated by the 
zeros of a certain function d{k) whose definition involves the initial conditions as well 
as the boundary values at x = 0. Thus, the boundary can have a nontrivial effect on 
the asymptotics of the solution. 

1.1 Purpose of the paper 

The aspects of the half-line spectral theory mentioned above have been analyzed in 
general terms in several papers (see for example [31 HJ El El HOl [El [THl [20l EH [22] ) , 
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but very few explicit examples have been studied either analytically or numerically. 
This is partly due to the fact that the effective construction of the solutions of IBV 
problems involves a certain step which is nonlinear: Despite the fact that the solution 
u{x, t) can be recovered from the solution of a RH problem formulated uniquely in 
terms of the values of u{x, 0) on the half-line {x > 0, t = 0} together with the 
values of u{x,t) and its x-derivatives on the boundary {x = 0,t > 0}, this does not 
provide directly an effective formula for solving the problem. Indeed, for a given IBV 
problem only part of the boundary values are prescribed as boundary conditions. The 
remaining boundary values cannot be independently specified, but are determined by 
the analysis of the so-called global relation. This yields the unknown boundary values 
in terms of the known initial and boundary conditions via the solution of a nonlinear 
Volterra integral equation [21 |TT1 [23| . 

Nevertheless, there exists a class of IBV problems for which the global relation 
can be solved explicitly. In this case all the ingredients needed for the formulation 
of the RH problem can be obtained as effectively as in the case of the initial value 
problem. These so-called linearizable IBV problems provide a class of problems for 
which the asymptotic behavior of the solution on the half-line can be analyzed with 
the same effectiveness as the corresponding problem on the infinite line. 

The purpose of the present paper is to provide the analysis of a number of concrete 
examples of linearizable IBV problems. The explicit formulas presented here elucidate 
the influence of the boundary on the asymptotic behavior of the solution. 

1.2 The KdV with dominant surface tension 

Our approach for constructing explicit solutions will be illustrated using the following 
equation: 

Ut + Ux-Uxxx + QuUx = 0. (1.1) 

This equation is the KdV equation in the case of dominant surface tension. Indeed, 
the Korteweg-de Vries equation for long one-dimensional, small-amplitude, gravity 
waves propagating in shallow water is 

where rj is the surface elevation of the wave above the equilibrium level h, a is a 
parameter related to the uniform motion of the liquid, g is the gravitational constant, 
T is the surface tension, and p is the density cf. pp. In the presence of sufficiently 



large surface tension, the parameter a is negative, and equation (1.2) can be brought 



into the nondimensional form (1.1) 



The reason for considering equation (1.1) instead of the standard KdV equation 



is that whereas the IBV for the standard KdV equation on the half-line is well-posed 



with one boundary condition at x = 0, equation (1.1) is well-posed if two boundary 
conditions are specified at x = 0. For our purposes the latter situation is favorable 
since it gives more freedom in the choice of boundary data at x = 0. In particular, 
equation (1.1) admits the following two classes of linearizable boundary conditions 
(see section |3 ) : 



4 



(a) n(0,t) =u^^(0,t) = 0, t>0, 

(b) n(0, t) = X, u^xiO, t) = x + ^X^, t>0, for some X S M with x 7^ 0. 

We will refer to boundary conditions of these two types as linearizable of type (a) or 
type (b), respectively. 



1.3 The examples 

We will consider three classes of examples: (i) Solitons (ii) IBV problems with ex- 
ponential intitial profiles (iii) IBV problems with box-shaped initial profiles. These 
classes of examples allow us to analyze problems with both types of linearizable 
boundary conditions: the IBV problems with box-shaped intitial profiles are associ- 
ated with linearizable boundary conditions of type (a), and the IBV problems with 
exponential intitial profiles are associated with linearizable boundary conditions of 
type (b)g 

We believe the study of explicit examples is of interest since: (a) it provides a non- 
trivial verification of the general formalism; (b) it reveals interesting properties that 
were previously unnoticed; (c) it makes more accessible the general rather abstract 
theory. 



1.3.1 Solitons 



Equation (1.1) formulated on the line, —oc < x < 00, admits solitons. We can 
construct an explicit solution to an IBV problem on the half-line by simply restricting 
such a soliton solution to x > and choosing compatible boundary conditions at 
X = 0. We will consider the restrictions of one-solitons and two-solitons and also 
comment on the case of rational solitons. 

The IBV problems constructed in this manner are of course artificial, but they 
are important for pedagogical considerations: In this case the spectral functions and 
all the ingredients of the RH problem can be computed directly, thus the analysis of 
these problems provides an important check on the general formalism and also gives 
us insight into nontrivial examples. 

We will compute the function d{k) whose zeros generate the solitons and verify 
that the general theory predicts the correct answer. Moreover, it turns out that for 
certain choices of the parameter values, these restricted soliton solutions satisfy a 
linearizable boundary condition at x = 0. Hence, for these cases we can also apply 
the approach available for linearizable boundary conditions. In this approach the 
function d{k) is replaced by another function A(fc) constructed in terms of the initial 
data alone; the theory claims that the zeros of A (A;) generate the asymptotic solitons. 
We again verify that the linearizable theory gives the correct result. 



^For the soliton solutions we can solve everything explicitly and will not need the special machinery 
available for linearizable boundary conditions. Nonetheless, for some choices of the parameters the 
boundary conditions turn out to be linearizable of type (b), and in these cases we will implement the 
linearizable approach as a check on the formalism. 
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1.3.2 Exponential intitial profiles 

Consider the family of initial profiles Er{x) defined by 

Er{x) = ^e"™ (r^ - l) , r > 0, r / 1, r constant. (1.3) 
The initial and boundary conditions 



Ur{x,0) = Er{x), X > 0, 

UriO, t) = Xr, UrxxiO, t) = Xr + ^Xr^ t > 0, Xr = ^ {r^ " l) , 



(1.4) 



define an IBV problem for equation (1.1) which is linearizable of type (b). 

In the case of linearizable boundary conditions, the formulation of the relevant RH 
problem requires only the analysis of the x-part of the Lax pair with the potential 
given in terms of the initial conditions. In general, this analysis involves solving 
a linear Volterra integral equation. However, for the present example, due to the 
simple form of the initial profiles Er{x), the eigenfunctions of the x-part of the Lax 
pair can be expressed explicitly in terms of hypergeometric functions. This leads 
(via the formalism for linearizable IBV problems) to an explicit expression for the 
function A(A;) (whose zeros correspond to the solitons present asymptotically in the 
solution Ur{x,t)). Even though the formula for A{k) is complicated, we can find the 
zeros numerically, and hence predict the number and shapes of solitons generated 



asymptotically by the initial and boundary data in (1.4). For example, we will find 
that for large r there exists asymptotically only one stationary (time-independent) 
soliton, while as r ^ the number of generated solitons appears to grow indefinitely. It 
is explained in section [5] that many of the relevant qualitative results can be explained 
in terms of the variations with respect to r of the shape of the initial profile and of 
the boundary conditions. 

1.3.3 Box-shaped intitial profiles 

Another natural and interesting class of initial profiles for which the x-part of the 
Lax pair can be solved explicitly is the class of box-shaped profiles defined as follows. 
Given data / = {h, L, xq} C M x M+ x M, we let 



Piix) 



0, < X < xo, 

h, Xq < X < Xo + L, 

0, Xo + L < X. 



The function /3/(x) represents a box of height h and length L positioned with its left 
side at X = Xo- The initial and boundary conditions 

fn/(x,0) = ^/(x), x>0, 
\uiiO,t) = ui^^iO,t) = 0, t>0. 



determine an IBV problem for equation (1.1) which is linearizable of type (a). In 



this case we can again derive an explicit expression for the function A(A;) and study 
numerically its set of zeros responsible for asymptotic solitons. 
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Note that the functions f3i{x) are well-defined also on the line and have decay as 
X — )• ±00 (this is in contrast to the case of the exponential initial profiles defined in 
(1.3), which do not decay as x — )• —00). Therefore, we may also apply the ordinary 



1ST to determine the asymptotic behavior of the solution of (1.1) with initial profile 
f3j{x) in the absence of a boundary. Naturally, we find that the number and shapes 
of asymptotic solitons in the absence of a boundary is independent of the position 
xq of the box. However, in the presence of a boundary we find that the number of 
solitons does depend on xq. The IBV problem with initial and boundary data given 



by (1.5) therefore provides a concrete example elucidating the effect of the boundary 



on the asymptotics. 



1.4 Outline 

Before analyzing particular examples, we review in section [2] the general theory of 
the 1ST for equation (1.1) on the line and on the half-line. Our aim is to make the 



presentation as unified as possible, so that the similarities (and differences) of the 
line and half-line problems become apparent. 

In section [3] we consider the special methods available for linearizable IBV prob- 
lems. The steps required for this analysis were already outlined in [12]; we will state 
the final results but also consider certain aspects in greater detail. 

In sections [4] - [6] we consider the three main classes of examples: Solitons; IBV 
problems with exponential intitial profiles; and IBV problems with box-shaped initial 
profiles. 

In the appendix we give the formulas necessary relating the matrix and scalar Lax 



pairs associated with equation (1.1). 



2 Spectral theory 



In this section we review the spectral theory of equation ( 1.1 ) on the line and on the 



half-line. A general introduction to the Inverse Scattering Transform can be found in 
PP for the case on the line, while the tools needed for the case on the half-line were first 
introduced in [51 [9] (see also [12] ) . Our treatment here seeks to implement the spectral 
theory for the line and half-line problems in a way that emphasizes the similarities 
between the two approaches. We will define four eigenfunctions of the Lax pair 



associated to (1.1); the eigenfunctions fJ.i,fJ,2,fJ-3 will be used in the analysis of the 



half-line problem, whereas the infinite line problem utilizes the eigenfunctions /is and 
/i4. For both the infinite line and the half- line problem, we then express the solution 



u{x,t) of equation (1.1) in terms of the solution of a 2 x 2 matrix Riemann-Hilbert 
problem. 



2.1 Lax pair 

Let 



1 

1 



0-2 



-i 

1 



0-3 



1 
-1 



7 



Equation (1.1) admits the Lax pair formulation 



I Hcc-ik[a3,fj] = Vin, ^^^^ 
\ ^j.t + i{k + 4k^)[a3, fi] = V2H, 

where e C is a spectral parameter, fi{x, t, fc) is a 2 x 2 matrix valued eigenfunction, 
and 

u 

Vi{x,t,k) = —{a2-ia3), (2.2) 
V2{x,t,k) = -2kua2 + UxCTi + ~^ — ^(^0-3 -0-2). (2.3) 



We can write the Lax pair (2.1) in differential form as 

d (^e'^-'''-'+^k+Ak^)t)-^^{x, t, A;)) = W{x, t, k), (2.4) 

where 

and (T3 acts on a 2 x 2 matrix A by a^A = [0-3, A],e"'^A = e'^^Ae~'^^. 
2.2 Bounded and analytic eigenfunctions 



Let equation (2.4) be valid for < t < T and —00 <x<oo(0<x<oo when 
considering the half-line problem) where T < 00 is a given positive constant. Unless 
otherwise specified, we suppose that the solutions are defined for —00 < x < 00; 
we will indicate what changes are necessary in the case of the half-line problem. 
Assuming that the function u{x, t) has sufficient smoothness and decay, we introduce 



four solutions /ij, j = 1,2,3,4, of (2.4) by 



Ax,t) 

fij{x,t,k) = I+ / e-'^-^^+^''+^'''^'^^''W{x',t',k), (2.5) 

where (xi,ti) = (0,T), (^2,^2) = (0,0), (x3,t3) = {00, t), and {x4,U) = (— oo,i). The 
function is only defined when considering the infinite line problem. If T = 00, the 
function fii is only defined if n(0,t) decays to zero as t — )■ 00. Since the one-form 



W is exact, the integral on the right-hand side of (2.5) is independent of the path 

of integration. We choose the particular contours shown in Figure |2.1[ This choice 
implies the following inequalities on the contours, 

(xi, ti) ^ {x,t) : x' - X <0, t' -t>0, 

(^2,^2) — 7- (x, t) : x' — X < 0, t' — t < 0, 
ix3,h) (x,t) : x' - X > 0, 
(x4, ti) — )• (x, t) : x' — X < 0. 
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(0,T) 


(x,t) 


>-(x,t) 














(x,t) 


(x,t) 



JX3 



1^4 



Figure 2.1 The contours of integration for the solutions fii, fi2, /U3, and m of (2.1) 



The second column of the matrix equation (2.5) involves exp[2i(— A;(x' — x) + {k + 
4/c^)(i' — t))]. Using the above inequalities it follows that this exponential is bounded 
in the following regions of the complex fc-plane, 



{xi,ti) - 


^ {x,t) 


{Imk > 0} n {Im 


[k + 4A;3] 


>0}, 


{X2,t2) - 




{ImA; > 0} n {Im 


[A; + 4A;3] 


<0}, 


(3^3,^3) - 




{ImA; < 0}, 






(3^4,^4) - 




{ImA; > 0}. 







We define the sets (see Figure |2.6| 

Di = {ke C\lmk < and Im [A; + 4A;^] > 0}, 
D2 = {ke C|Im A; < and Im [k + 4A;^] < 0}, 
Ds = {k£ C|Im A; > and Im [k + Ak^] > 0}, 
L>4 = {A; G C|ImA: > and Im [A; + 4A;^] < 0}. 

The second column vectors of fii, /X2, /U3, /X4 are bounded and analytic for A; € C 
such that k belongs to -D3, -D4, Di U D2, and U 1)4, respectively. We will denote 
these vectors with superscripts (3), (4), (12), and (34) to indicate these boundedness 
properties. Similar conditions are valid for the first column vectors. We obtain 



(2) (3) 



(1) (4) 



^J'3 



(34) (12) 
A''3 ' /^3 



fj,4 



(12) (34) 
A''4 ' A*4 



We deduce from (2.5) that /ii is an entire function of k, and that /i2 is an entire 



function of A; if T < 00. 
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Let us collect some properties of the eigenfunctions of (2.1 ). As A; — )• 0, a solution 



^{x,t,k) of (2.1) satisfies 
li{x, t, k) 



ia{x,t) / 1 1 



k 



-1 -1 



+ 0(1), k^O, 



(2.6) 



where a{x, t) is a real- valued function cf. jl2j . From the form of the matrices Vi and 
V2 defined in ( |2.2[ ), it follows that each eigenfunction fj, of (2.1) obeys the following 
symmetries: 



Hii{x, t,k) = ^22{x, t,k), 
fiii{x,t,k) = H22ix, t,-k) 

It also holds that 

Hj{x, t, k) = I + O 
Furthermore, note that 



fii2{x,t,k) = fi2i{x,t,k), 
IJ.i2{x,t,k) = H2iix, t,-k). 



/c 00, j = 1,2,3,4. 



(2.7) 
(2.8) 



det^j = l, i = 1,2,3,4, 



(2.9) 



and that any two solutions /i and /i of (2.4) are related by an equation of the form 

fi{x, t, k) = /i(x, t, fc)e-^(-'=^+('=+^'=')*)^3Co(A:), (2.10) 

where CQ{k) is a 2 x 2 matrix independent of x and t. 

The /ij's are the fundamental eigenfunctions needed for the formulation of a 
Riemann-Hilbert problem in the complex /c-plane. The eigenfunctions /xi,/i2,^3 are 
used for the half- line problem, while H3,fJ-4 are used for the infinite line problem. 

2.3 Spectral theory on the Hne 



We now restrict attention to the infinite line problem. In view of (2.10), we may 
define a spectral function s^{k) bjj^ 



fi3{x,t,k) = ^4(x,t,A:)e-*(~^^+(^-+^'=')*)^^s^(A;), 



Evaluation of (2.11 ) at t = and x — )• — 00 gives 



s^{k) = I 



(yiH3){x,0,k)dx, 



Im k = 0. 



Imk = 0. 



(2.11) 



(2.12) 



Moreover, by (2.9), 



del s^{k) = 1. 



We infer from the symmetries (2.7) that there exist functions a^{k) and b^{k) such 

that 

V(^) b^{k)\ 



s^{k) 



b^{k) a^{k) 



(2.13) 



■^Some quantities related to tlie problem on the line will be denoted with the superscript L to 
distinguish them from their half-line counterparts. 



10 



From the explicit expression (2.12 ) for s^{k) and the fact that the second column of /^a 
is defined and analytic in D1LID2, we deduce that a^{k) has an analytic continuation 



to Di U D2. It follows from (2.6) and (2.11) that 



+ 0(1), b'^ik) 



+ 0(1), 



A; ^ 0, 



for some real constant a. 



2.3.1 Residue conditions 

We assume that a^{k) has N simple zeros {fcj}^^ in the lower half-plane. These zeros 
automatically lie on the imaginary axis as we consider only real-valued solutions of 
(1.1). This can be seen as a consequence of the fact that the operator defining the 
x-part of the scalar Lax pair in (A.l) with eigenvalue fe^ is self-adjoint and so has a 
purely real spectrum. 



The second column of equation (2.11) is 



Im k = Q. 



(2.14) 



Applying det ( ■ ) to this equation and recalling (2.9), we find 



dei[ji^l'^\x,t,k),ijL\""'\x,t,k)'j =a''{k), /c € -Di U D2, 

where we have used that both sides are well-defined and analytic in the lower half- 
plane to extend the above relation to D1UD2. Hence, if a^{kj) = 0, then n^^'^\x, t, kj) 
(12) 

and //g (x, t, kj) are linearly dependent vectors for each x and t. It follows that there 
exist constants bj such that 

fi'i^\x,t,kj) = bje^'^-''^-^+^''^+^''f^'^fi^i^\x,t,kj), xeR,t>0. (2.15) 



,(12) 



Recalling the symmetries (2.7), the complex conjugate of (2.15) is 
Consequently, the residues of fi^^"^^ /a^ and fj,^^^^ /a^ at kj and kj are 



Res 



(12). (12). , X 

^4 l^,'^,ft^J _ ^4 l^,t,/CjJ _ 2j(-fcjX+(fcj+4fcf)t), (12)/ , N 



a^ik) 



Rcs ^^ ''(^,^,^) _ '{x,t,kj) _ (^.g-2j 



,(34), 



h a^{k) 



d^{kj 



i-kiX+(ki+4k^)t) (34)/ , r n 
^ ' ^ ^' '[X,t,kj), 



where a 



^ - '^"'^ and C 



dk 



3 - a^kj)- 
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Remark 2.2 In the traditional approach to the Inverse Scattering Transform the x- 



part of the Lax pair ( 2.1 ) is first analyzed at each fixed time t to yield time-dependent 



scattering data consisting of two functions a^{k, t) and b^{k, t), together with the ze- 
ros {kj{t)}i of a^{k, t) and the corresponding normalization constants {Cj{t)}^ . The 
t-part of the Lax pair is then used to determine the time-evolution of the scattering 
data. In our present formulation the x- and t-parts of the Lax pair are analyzed si- 
multaneously, which implies that no separate analysis is necessary in order to find the 
time-evolution. The relation between the scattering data in the traditional analysis 
{a^{k,t),b^{k,t),{kj{t)}^,{Cj{t)}^} and the data {a^ik),b^{k),{kj}^ ,{Cj}^} in 
our current treatment is 



a^ik,t) = a^{k), b^{k,t) = b^{k)e 



2i{kj 



-4kpt 



kj{t) = kj, Cj{t) 

In our current treatment the time-dependent exponential factors g='=2*(fcj+4fc^)* g^^g g^^, 
tomatically built into the equations (2.11) and (2.15). We have chosen the present ap- 
proach for two reasons: (i) It is the proper approach for dealing with initial-boundary 
value problems, so that adopting it also for the problem on the line makes the link 
to the half-line problem more visible, (ii) We believe it is more natural also for the 
problem on the line. 



2.3.2 Riemann-Hilbert problem 



We now describe how the solution to the initial- value problem for equation (1.1) 



on the line can be expressed in terms of the solution M{x, t,k) of a 2 x 2 matrix 
Riemann-Hilbert problem. 

By substituting the expansion 



777,(1) j^(2) / 1 



k — )• oo, 



in the x-part of (2.1), we find by considering the terms of 0{l/k) that 

u{x) 



-2im 



(1) 

22a;- 



Algebraic manipulations show that relation (2.11 ) can be rewritten in the form of the 
RH problem 

M^{x,t,k) = M+{x,t,k)J{x,t,k), Imk = 0, 
where the matrices M_, M^, J are defined by 



' (12) ^ 
„{12) 



(34) \ 
(34) f^i \ 



a^ik) ) 

1 



Imfc < 0: 



Im > 0: 



(2.16) 



J 



b^{k) 2i(-kx+(k+4:k^)t) 



b^jk) ^-2i(-kx+{k+ik^)t)^ 
a^{k) 

1 



lm/c = 0. (2.17) 



a'-{k)aL(k) 
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The contour for this RH problem is the real axis. 

We summarize our discussion of spectral theory on the line in the following theo- 
rem, which reduces the Cauchy problem for equation ( 1.1 ) to a 2 x 2 matrix Riemann- 
Hilbert problem. We use the notation [A]i ([^]2) for the first (second) column of a 
2x2 matrix A. 



Theorem 2.3 Given initial data uo{x) such that uq{x),xuq{x) G L^(M"'"), let Vi 
j^{o'2 — ics) CLnd let fi3{x, 0, k) and fi4{x, 0, k) be the unique solutions of the followi 
Volterra linear integral equations: 

/•QO 

^is{x,0,k) = I- / e-"'^-'-^^'^^Vifi3){x',0,k)dx', 

J X 

^li{x,0,k) = I + e-*^(^'-^)^3(Vi/X4)(x',0,/fc)dx'. 
Define {a^{k),b^{k),Cj} by 

[s^ik)]2, s^{k) = I- i e-'^''^-'{Vifi3)ix,0,k)dx, lmk = 0, 



oo 



a^{k) 



oo 



(2.18) 



and 

[/i4(x, 0, kj)]i = d^{kj)Cje-^'^^''[fi3ix, 0, kj)]2, j = l,...,N, (2.19) 

where we assume that a^{k) has N simple zeros {kj}jLi C iM^. 
Then 

• a^{k) is defined for Imk < and analytic in Imk < 0. 

• b^{k) is defined for /c G M. 



• a^{k)aL{k) - b^{k)b^{k) = 1, ke 

• a^{k) = a^{-k), Imk < 0. 



Moreover, the solution u{x,t) of equation (1.1) with initial data uq{x) is given by 

u{x, t) = -2i lim kd^ {M22{x, t, k)) , (2.20) 

k—^oo 

where M(x, t, k) is the unique solution of the following RH problem: 

M[x,t,t,)-<^ M_(x,t,fc), lmk>0, 
is a sectionally meromorphic function. 



M^{x,t, k) = M^{x,t,k)J{x,t,k) fork G M, where J is defined in (2.11). 
M{x, t, k) has the asymptotic behavior 

Mix,t,k) = I + (^j , k^oo. (2.21) 
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For some real function a{x,t), 



M+{x,t,k) 
M_{x,t,k) 

For some real constant a, 
ia 



ia{x,t) fo 1 



k \0 -Ij' 
ia{x,t) / I 0\ 

~^V-i or 



k^O 
k-^0. 



(2.22) 
(2.23) 



a^(fc) = ^ + 0(l), bHk) = ~ + Oil) 



0. 



The first column of Mj^ has simple poles at k = kj, j = 1,...,N, and the 
second column of Af_ has simple poles at k = kj , j = 1, . . . , N . The associated 
residues are given by 



Res[M{x,t,k)]i =Cje2*(-'=^-^+('=:'+4^|)*)[M(x,t,A:j)]2, 

kj 

Res[M{x,t,k)]2 =Cje-2^(-^^'^'+(^^+^^'W[M(x,t,^j-)]i, 



j = l,...,N, (2.24) 
j = l,...,N. 

(2.25) 



Remark 2.4 (Regularized RH problem) The RH problem of Theorem 2.3 has a 



singularity at A; = (see equations (2.22) and ( |2.23 )). It is sometimes more convenient 
to consider the regularized RH problem with no singularity at k = 0. The solutions of 
the two problems are related as follows: If M{x, t, k) is a solution of the regularized 
RH problem (i.e. M satisfies the jump and residue conditions, but is regular at 
= 0), then 



M{x,t,k) 



1 + 



iy{x,t) ( 1 



2k 



M{x,t,k), y{x,t) = 2i lim {kMi2{x,t,k)), 

k—^oo 

(2.26) 



is a solution of the RH problem of Theorem 2.3 cf. 



Remark 2.5 (Symmetries of the solution of the RH problem)) It follows from 



the symmetries (2.7) that the solution M{x,t, k) of the Riemann-Hilbert problem in 



Theorem 2.3 respects the symmetries 



Mii(x,t,A;) = M22(x,i,fc), M2i{x,t,k) = Mi2{x,t,k), (2.27) 
Mii{x,t,k) = M22{x,t,-k), Mi2{x,t,k) = M2i{x,t,-k). 



Moreover, if one enforces (2.27), then only one of the two residue conditions (2.45)- 



(2.46) needs to be verified since the other condition is a consequence of symmetry. 
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2.3.3 Soliton solutions 



We next outline how to find the solitons of equation (1.1) and the corresponding 



eigenfunctions of the Lax pair (2.1). We will follow these steps in subsections 4.1 and 



4.2 below to compute the one and two-solitons and their associated eigenfunctions. 

The solitons correspond to spectral data {a^{k),b^{k),Cj} for which b^{k) van- 
ishes identically. In this case the jump matrix J in ( 2.17| ) is the identity matrix and 
the RH problem of Theorem 2.3 consists of finding a meromorphic function M{x, t, k) 
satisfying ( |2.21[ ), the residue conditions ( |2.24[ )-( |2.25] ), and the conditions ( |2.22[ )-( |2.23| ) 
at k = 0. We first construct the solution M to the regularized RH problem with no 
singularity at = 0; subsequently, the solution M is found from (2.26). In what 



follows we will suppress the (x, i)-dependence of M and write M{k) for M{x, t, k). 



Let Oj = -kjx + {kj + 4k])t. From (2.21) and (2.24) we find 



[Mik)] 



+ 



EG,e 



2iei 



[Mik 



(2.28) 



If we impose the symmetries (|2.27|), equation (2.28) can be written as 

.^Cj^fMuik 
,0y ^ ^^k-kj \M22{kj)^ 



( M22{k) \ 
\Mi2{k)) 



(2.29) 



Evaluation at A;„ yields 



,Mi2ik„ 



+ 



h 



_ (Mi2{kj) 
kj \M22{kj) 



n 



.N. 



(2.30) 



Solving this algebraic system for Mi2{kj) and M22ikj), j = 1, . . . , N, and substituting 
the result back into (2.29) gives an explicit expression for [M(A;)]2. The first column 



of M{k) is obtained by symmetry. This construction provides the solution of the 
regularized Riemann-Hilbert problem. The solution M of the original RH problem is 



found from (2.26), and the soliton solution u{x,t) is obtained from (2.20). 



2.4 Spectral theory on the half-hne 

We now turn to the spectral analysis for the half-line problem. The presentation will 
be brief; see for further details. We define s{k) and S{k) by the relations 

/i3(x, t, k) = /X2(x, t, /fe)e-^(-'=^+('=+^'=')*)^3s(fc), (2.31) 
/ii(x, t, k) = /X2(x, t, /fe)e-^(-'=^+('=+^'=')*)^-^5(A;). (2.32) 

Evaluation of (2.31) and ( 2.32[ ) at (x,t) = (0,0) and (x,t) = (0,r) gives the expres- 



sions 



s{k) = ;U3(0, 0, k), S{k) = Aii(0, 0, k) = f e*(^+^'=')^^V2(0, T, k) 



(2.33) 
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where the final equahty is valid only if T < oo. We use the following notation for s 
and S: 

/a(^ b{k)\ s(k)=(^ (2 341 

^ ^ \b{k) aik) J' ^ ' \B{k) A{k) J' ^ ' ^ 

The spectral functions a{k) and b{k) have the following properties: 

(i) a{k) and b{k) are continuous and bounded for k G DiU D2 and analytic in the 
interior of this set. 

{ii) a{k) = 1 + 0(1 /k), b{k) = 0{l/k), k ^ 00, /c G A U L>2- 

(iii) |a(A:)|2 - \b{k)\'^ = 1, A; G U D2. 

(iv) For some real constant a, 

trv try 

a{k) = ^ + 0{l), 6(fc) = -^ + 0(l), k^O. (2.35) 

The spectral functions A{k) and B{k) have the following properties: 

(i) A{k) and B{k) are entire functions bounded for k G Di U -D3 when T < 00. If 
T = 00, the functions A{k) and B(k) are defined only for k G Di U D^. 

{ii) A{k) = l + 0{l/k), B{k) = 0{l/k), k^oo, keDiUD^. 

c, r < 00, 

D^ U As, T = 00. 



(iii) ^(A;)^(A;) - B{k)B{k) = 1, G 



(iv) For some real constant /3, 

^(A;) = ^+0(1), B{k) = --^+0{l), k^O. (2.36) 
2.4.1 Global relation 

Applying Stokes' theorem to the domain {0 < x < 00, < t < T} and the closed 
one- form W with /i = one finds the following so-called global relation (see |12j): 

[0 for keDi, T = oo, 

where 

/>oo 

c+(A;)= / e-^''-'{V^^,s)i2{x',T,k)dx'. 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 2.6 The contour for the Riemann-Hilbert problem in the complex k-plane. 



2.4.2 Riemann-Hilbert problem 



Equations (2.31) and (2.32) can be rewritten in the following form, expressing the 
jump condition of a 2 x 2 RH problem: 

M_{x,t,k) = M+{x,t,k)J{x,t,k), keDiDDj, i,j = l,...,A, 

where the matrices M_, M-|_, and J are defined by 



M4 



M4 



, (12)' 

,a{ky^-^ ^ 



d{k), 



, k€Di; 
, k£D3; 



M_ 



M_ 



,d{kr^^ 



a{k)_ 



k G D2; 
k £ D4; 



d{k) = a{k)A{k) -h{k)B{k), k £ D2; 



J{x, t, k) 



keDinD2 



J2 = J3Ji Jl k£D2nD3 



with 



k e DsHDi 
k e DiDDi 



a(k) 

1 

a(fc)a(fc) 



(2.38) 



(2.39) 
(2.40) 



1 -r(fc)e-2*^('^) 

1 

(2.41) 
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e{x,t,k) = -kx+ {k + 4:k^)t; T{k) 



B{k) 



a{k)d{ky 

The contour for this RH problem is depicted in Figure [2^ 



keD2 



(2.42) 



such that kj G Di, j = l,...,ni, and 



Assumption 2.7 We assume that 

(i) a{k) has simple zeros {kj}jLi C il 
kj £ D2, j = ni + l,...,N. 

(ii) d{k) has A simple zeros C iM~, such that Xj G D2, j = 1, . . . , A. 

(iii) None of the zeros of a{k) coincides with a zero of d{k). 

Theorem 2.8 Letuo{x), x > 0, and {gj{t)}i, t G [0,T), be such that uq{x) , xuq{x) G 
L^(M+) and gj{t),tgj{t) G L^{[0,T)), j = 1,2,3. Define the spectral functions a{k), 
h{k), A{k), and B{k) according to (2.33) and (2.34), where fii{0,t,k) and fi3{x,0, k) 
are obtained as the unique solutions of the Volterra linear integral equations 

rt 



fii{0, t,k) = 1+ [ e 
Jt 

/•oo 

H3{x,0,k) = I 



^H-'-)^^^Vifis){x',0,k)dx', 



(2.43) 
(2.44) 



and Vi{x, 0, (), V2(0, t, Q are given by equation (2.2) in terms of the initial and bound- 
ary values 

no(a;) = u(x,0), goit) = u{0,t), gi{t) = UxiO,t), g2{t) = Uxx{0,t). 
Suppose that the initial and boundary values are compatible in the sense that 



• they are compatible with equation (1.1) at x = t = 0. 

• the spectral functions satisfy the global relation ( '^.gTp p] 

Assume that the possible zeros {kj}^ of a{k) and {Xj}i of d{k) are as in Assumption 



2.1 . Define M{x,t, k) as the solution of the following 2x2 matrix RH problem: 



M is sectionally meromorphic away from the boundaries of the Dj 's, j 



1,...,4. 



M satisfies the jump condition 

M^{x,t,k) = M+{x,t,k)J{x,t,k), keD.DDj, i,j = l,...,4, 

where M is M_ for k £ D2U D^, M is M_|_ for k € Di U D3, and J is defined 
in terms ofa,b,A, and B by equations (2.39)-(2.42). 



If T < 00 we require that the function c on the right-hand side of (2.371 be continuous and 
bounded for ( £ DiU D2, analytic in Di U D2, and c^(fc) = 0{l/k) as A:—)- cxd. 
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The first column of M has simple poles at k = kj, j = 1, . . . , ni, and k = Xj, 
j = 1, . . . , A. The second column of M has simple poles at k = kj and k = \j, 
j = 1, . . . , A. The associated residues satisfy the following relations: 

Res[M(x,t,fc)]i = I e^'^fe)[M(x,t,fc,)]2, j = l,...,2ni, (2.45) 
kj a[kjjb(kj) 

Res[M{x,t,k)]2= ^ e-"^^fe)[M(x,t,fe,)]i, j = l,...,2ni, (2.46) 

Res[M{x,t,k)]i =Resr{k) e'^'^^^^^[M{x,t, Xj)]2, j = 1,...,2A, (2.47) 



Res[M(x,t,A:)]2 =Resr{k) e-^'^^^^\M{x,t, Xj)]i, j = 1,...,2A, (2.48) 



where 9{kj) = —kjX + {kj + 4:kj)t, and Res-^. T{k) denotes the residue of the 
function k i— )■ T[k) at k = Xj. 

• M{x, t,k) = I + (l) , k^ CO. 

• For some real function a{x,t), 

M(x,t,fc)~^^(^° _\), k^o,keD,. 

Then M{x, t, k) exists and is unique. 

Define u{x,t) in terms of M{x,t, k) by 

u{x, t) = -2i lim kdx {M22{x, t, k)) . (2.49) 
Then u{x, t) solves equation Furthermore, 

u{x,Qi) = uq{x), u{0,t) = go{t), UxiO,t) = gi{t), and u^xi^^t) = g2{t). 



Let us also point out that Remark 2.4 applies also to the RH problem for the 



half-line problem: In order to find the solution M, we may first find the solution M 
of the regularized RH problem with no singularity at A; = 0, and then recover M from 



equation (2.26) 



3 Linearizable boundary conditions 



It was shown in Theorem |2.8| that the solution u{x,t) of equation (1.1) on the half- 
line can be expressed through the solution of a 2 x 2 matrix RH problem, which is 
uniquely formulated in terms of the spectral functions a{k), b{k), A{k), and B{k). 
The functions a{k) and b{k) are defined in terms of the initial data uq{x) through 



the solution of the linear Volterra integral equation (2.44). However, the spectral 



functions A{k) and B{k) are, in general, not as readily obtained: The construction 



of A{k) and B{k) via the linear Volterra integral equation (2.43) requires knowledge 
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of go{t), gi{t), and 52 (i), whereas the given boundary conditions impose only two 
conditions among these three functions; the additional condition needed to determine 



go{t), gi{t), and g2{t) is the requirement that they satisfy the global relation (2.37). 
In general, this problem involves solving a nonlinear Volterra integral equation. 

However, for a particular class of boundary value problems it is possible, using 
only the algebraic manipulation of the global relation, to compute functions A{k) and 
B(k) which effectively replace A(k) and B(k). More precisely, the solution M to 
the RH problem associated with A, B instead of A, B, can be directly related to the 



solution of the original RH problem (see Theorem 3.1 below). In particular, u{x,t) 
can be recovered from the large k asymptotics of M. When T = 00 the functions A 
and B coincide with A and B, so that the two RH problems are identical. The class of 
boundary value problems which yield to this approach are referred to as linearizable. 
Thus, for linearizable IBV problems the half-line formalism is as effective as the 
formalism on the line. An analysis of linearizable boundary conditions for equation 



( 1.1 ) was outlined in [12j . However, our consideration of specific examples will require 
some additional information in order to keep track of specific branches of solutions 
to third-order polynomials; we will therefore in this section present a more detailed 
version of the analysis of [12j. 

The derivation involves considering for each k the roots j^(/c) of the polynomial 

iy{k) + Aiy{kf = k + 4k^. (3.1) 

For each k, there are three roots I'jik) £ C, j = 1,2,3 (one of these roots equals 
k), but it is impossible to choose a consistent numbering of these roots on the whole 
complex plane. Indeed, the two roots not equal to k get interchanged as k encircles 
one of the points k = ibi/2. However, it is possible to fix a numbering of the roots for 
k lying in the restricted set Di U . For /c G Z)i U D3 one root lies in Di while the two 
other roots lie in D^. We may thus define continuous functions z^i : Di U 1)3 — )• Di 
and Uj : Z)i U 1)3 — )• -D3, J = 2, 3, such that 

Uj{k) + 4i^j{kf = k + 4k^, keDiUD^, j = 1,2,3. 

For each k G Di U exactly one of the complex numbers Vjik), j = 1, 2, 3, equals k. 
We extend the domain of definition of Vj, j = 1, 2, 3, to include D2 U D4 by defining 



Vj{k) = Uj{k), k(^D2UDi. 
With this notation, the basic result is the following. 



Theorem 3.1 Let u{x,t) satisfy equation (1.1) on {x > 0, < t < T} with T < 00 
together with the initial condition 

u{x, 0) = uq{x), < X < 00, 

and one of the following boundary conditions: 

(a) n(0,t) =u^^(0,t) = 0, 
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(b)u{o,t) = x, u^xio,t) = x + ^x^ xe 

In the case (a) we define 



A{k) =a{k)±M-b{k), 



k G D2, 



biulik)) 



a{k)A{k)' 
In the case (h) we define 

vi{k) + k 



kGDo 



h{k) 



ui{k) — k 



Akvi{k)\ 
X J 



A{k) = a{k) 



h{k)a{ui{k))-b{ui{k)) 
h{k)b{vi{k)) - a{u^{k)) 

' keD2. 



b{k), k£D2 



a{k)A{k)' 



(3.2) 
(3.3) 

(3.4) 
(3.5) 
(3.6) 



Assume that the initial and boundary conditions are compatible at {x,t) = (0,0). // 
T = 00, we assume that the boundary conditions are of type (a). Furthermore, assume 
that the following are valid: 



(i) a{k) has N simple zeros {kj}jL-^ C iM such that kj € Di, j = 1, 
kj £ D2, j = ni + l,...,N. 

(a) A{k) has A simple zeros {Aj}^ C iM^ such that Xj G D2, j = 1, . . . , A. 

(Hi) None of the zeros of a{k) coincides with a zero of A{k). 



,ni, anc 



Then the solution u{x,t) is given by equation (2.49) with M re placed by M , where 
M is the solution of the Riemann-Hilbert problem in Theorem 2.8 with jump matrices 
and residue conditions defined by replacing T{k) in (2.4I) and (2.4T)-(2.48) with 

Let Ji, i = 1,...,4, be the jump matrices defined according to {2.41) with T 
replaced by T. Then M is related to the solution M of the original RH problem by 

Mi = Mi, M2 = M2Ji^Ji, Ah = M^J^J^^, Mi = M4- (3.7) 

Moreover, if T = 00, 



r{k) = T{k), 

so that the two RH problems coincide. 



{3.i 



Proof. Suppose that it is possible to compute explicitly in terms of the given boundary 
conditions, a nonsingular matrix Nj{k) such that, for j = 1, 2, 3 and k £ C, 

[i{k + 4k^)a3 - ¥2(0, t, Uj{k))] Nj{k) = Nj{k) [i{k + 4A;3)(T3 - ^2(0, t, k)] . (3.9) 
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If Nj{k) exists and T < oo, it follows from ([232]), and K9\ that 



k £ 



(3.10) 



It turns out that matrices Nj(k) satisfying (3.9) exist for the boundary conditions in 
(a) and (b). 

We can now explain the idea of the proof: The spectral functions A{k) and B{k) 



enter the formulation of the RH problem in Theorem 2.8 only through the function 



T{k) defined in (2.42). The goal is to use the algebraic relations between S{vj{k)) 



and S{k) obtained from (3.10) together with the global relation (2.37) to express T{k) 



in terms of a, b, and the matrices Nj. This turns out to be possible when T = oo. 
When T < oo we will define modified spectral functions A{k) and B{k) and define F 
by the same equation (2.42) as we used to define T, but with A and B replaced by A 
and B, i.e. 



m 



B{k) 



k £ D2. 



(3.11) 



aik){a{k)A{k) - b{k)B{k)) 
It can then be shown th at t he solutions of the RH problems formulated in terms of 
r or r are related as in (3.7). 



The first step is to show that r(A;) as defined in (3.11 ) can be expressed as in (3.3 ) 



and (3.6) when u satisfies boundary conditions of type (a) and type (b), respectively. 



We consider the two cases below. 



3.1 Boundary conditions of type (a) 



Assume first that T < oo. In this case (3.9) holds with Nj = I, j = 1,2,3, which 



using (3.10) leads to 

A{u,{k)) = A{k), B{v,{k)) = B{k), i = 1,2,3. 



In view of this symmetry, the global relation (2.37) yields 

B{k)a{vi{k)) - A{k)b{ui{k)) = e2^(^+4'=')^c+(z^i(A;)), k £ DiU 



Moreover, (3.12) together with the definition of d, 

d{k) = a(^A{k) -b(k)B{k), k e D^, 

imply 



(3.12) 
(3.13) 
(3.14) 



d{uj{k)) = a{uj{k))A{k) -b{uj{k))B{k), keDiUD^, j = 2,3. (3.15) 



Equations (3.13) and (3.15) can be solved for A and B to give, for j = 2,3 and 



k S Di U D3, the following expressions: 

a(i/i ( /fc) )d{^^) + e2'(^+4*^')'^c+ {ui {k))b(UJ{I^) 



A{k) 



B{k) 



A,ik) 



b{vi{k))d{uj{k)) + e^^^^+^^)^c+{ui{k))a{uj{k)) 



(3.16) 
(3.17) 
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where 



Aj{k) = a{i^i{k))a{i^j{k)) - b{vi{k))b{uj{k)), k e DiU D3. 



The last terms on the right-hand sides of equations (3.16) and (3.17) involve the 
unknown function c^. We are therefore led to consider the RH problem with A and 
B replaced by 



A{k) 



A,{k) 



h{vi{k))d{Uj{k)) 

A,{k) 



(3.18) 
(3.19) 



Let r be defined by (|3.11|), thus 

m 



B{k)/A{k) 



a{k) a{k) - b{k)B{k) / A{k) 



a{k)A{ky 



k e D2, 



where A(A;) is given by (3.2). 

For T = 00 the derivation is similar, but simpler. In p articular, A and B coincide 
with A and B in this case. This establishes formula (3.3) for T in case (a). 

3.2 Boundary conditions of type (b) 

Let 



Then (3.9) is satisfied for 

1 



J = 1,2,3. 



N,{k) 



1 m) ■ 



Equation (3.10) gives, for /c G C (we again assume T < 00), 
if,{kf - l)A(z.,(A;)) = f,ik) [Bik] + f,ik)A{k)] - e^^^^+^>^')'r 

{f,{kf - l)B{v,{k)) = f,{k) [A{k) + f,{k)B{k)] - 



A{k) + f,{k)B{k)\ , 
(3.20) 

(3.21) 



Note that fj{k) = fj{k), j = 1, 2, 3. By taking the Schwarz conjugate of (3.20)-(3.21 ) 
we find, for k £ C, 



{f,{kY - l)A{u,{k)) = f,{k) B{k) + f,{k)A{k) 



-2i{k+Ak^)T 



[A{k) + f,{k)B{k)] 



ifjiky - l)B{u,{k)) = f,{k) A{k) + f,{k)B{k) 



(3.22) 

-2i{k+ik^)T [B{k) + fj{k)A{k)] . 

(3.23) 
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Eliminating B{k) + fj(k)A(k) from equations (3.21) and (3.22), we find 



o2i 



+ T77T + ^(^)' ^^C' i = l'2,3. (3.24) 
fj[k) 



Similarly, eliminating A{k) + fj{k)B{k) from equations (3.20) and (3.23), we find 

e^i(k+ik^)T B{uj{k)) B{k) 



A{u,{k)) 



+ 



+ A{k), keC, j = 1,2,3. (3.25) 



Mk) ' Mk) 

In order to obtain A(k) and B{k) for A; G 1)3 we replace A{i/i(k)) and B{vi[k)) 
by the right-hand sides of (3.24) and (3.25) in the global relation 

B{ui{k))a{ui{k)) - A{vi{k))h{ui{k)) = e2^('=+4^-')*c+(i/i(A:)), k G D^. (3.26) 

This gives 



m 



a{vi{k)) 



h{k) 

g2i(fc+4A:3)T 



6(i.i(fe)) 



a(z^i(/c)) 



A(z/i(fc))a(z.i(A:)) - B{vi{k))b{ui{k)) 



Using the definition (2.39) of d, which is valid for k G D3, we infer that 

'KMk)) 



m 



a{vi{k)) 
Mk) 



HMk)) 



B{k) 



Mk) 



c+{Mk)) + 



d{vi{k)) 
Mk) 



a{vi{k)) 



k £ D^. 



Solving this equation together with (3.14) for A and B, we find 
1 



Mk) 



A,{k) 



aiMk))-'W^)dik) 



Mk) 

+ e'^^'^^'^'^^m { cHMk)) + 

Ji\k) 



(3.27) 



kGDs, 



B{k) 



Ai{k) 



HMk))-'^)dik) 



Mk) 



(3.28) 



+ e 



2i{k+Ak^)T 



a{k)\c+{ui{k)) + 



d{Mk)) 
Mk) 



k G ^3, 



where, for k £ D2, 



Ai{k) = a{k)a{vi{k)) - b{k)h{vi{k)) + 



Mk) 



-a{k)b{ui{k)) + h{k)a{vi{k)) 
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Some of the terms of the right-hand sides of equations (3.27) and (3.28) involve 



unknown functions, so that we are led to consider the RH problem for which A and 
B are replaced by the following expressions: 



A{k) 
B{k) 



1 



Then 



Ai(A:) 
1 

Ai(fc) 



a{ui{k)) 



fiik) . 
h{k) 



dik), 



d{k), 



k£D3, 
keDs. 



(3.29) 
(3.30) 



B{k)/A{k) 



a{k) a{k)-b{k)B{k)/A{k) 
^ keD2, 



a{k)A{k)' 



where A(A;) is given by (3.5). This establishes formula (3.6) for F in case (b). 



It remains to show that the solu tions M and M of the RH problems defined in 
terms of F a nd F are related by (3.7), and that T(k ) = T (k) when T = oo. It will then 



follow from (3.7) that u{x,t) is given by equation (2.49) with M replaced by M. The 
arguments necessary to prove (3.7) are presented in [9|. The equality F(A;) = T{k) 
when T = oo is a consequence of the fact that A = A and i? = ^ in this case. The 
restriction to case (a) when T = oo is necessary in order for the eigenfunction fii 



(and thus also for F) to be well-defined, see subsection 2.2 This completes the proof 
of Theorem 13.11 □ 

Remark 3.2 (Characterization of linearizable boundary conditions) The lin- 
earizable boundary conditions can be characterized as those boundary conditions 
which are time-independent and for which the total mass Jq°° udx is conserved. This 
follows directly from the fact that 



dt 



u{x, t)dx = \uxx — u — Su^] 



x=0' 



Equation (3.31) is a consequence of the conservation law 

Ut + {u- Uxx + 3ii^)^ 



0. 



(3.31) 



(3.32) 



4 Soliton profiles 



We now begin the study of particular examples of IBV problems. The first exam- 



ples are derived from soliton solutions of equation (1.1). The expressions for these 



solutions and the corresponding eigenfunctions of the Lax pair can be found from 
the RH problem in Theorem |2.3| as indicated in subsection 2.3.3 We will consider 



the one-soliton and the two-solitons, and will also comment on the case of rational 
solitons. 
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4.1 One-solitons 



In order to find the one-soliton of equation ( 1.1 ) and the corresponding eigenfunctions, 
we consider the algebraic system (2.30) with N = 2 and k2 = —ki (in view of the 
symmetries (2.7), if a^(fc) = then also a^{-k) = 0). Since 9j = -kjX + {kj + Ak^-)t, 

we have 9i = —^2- The identity a^{—k) = a^{k) implies that d^(/c2) = —a^{ki). 
Moreover, (2.15) together with the symmetries (2.7) imply that 62 = bi. Hence, 
C2 = —Ci and, using (2.27), the system (2.30) can be written as 



( M22{ki) 
\Mi2{ki) 



h - ki \M22{ki)J ki + h \M22ik1] 



(4.1) 



Since we consider only real- valued solutions, the zeros of will lie on the imaginary 
axis, so we may assume that ki = ip = k2 where p > 0. We parametrize the 
normalization constant Ci as Ci = —iape'^^°, where a = ±1 and xq G M. Fr om ( 2.27) 
we deduce that Mi2{ki) and M22(^i) are real. Substituting the solution of (4.1) into 
(2.29), we obtain 



Muik) 

M22{k) 



2ae 



2[itp^+{x-t)p+xo) 



-1 + e4(4tp3+(x-t)p+xo)^ (p _ ' 

2g4(4tp3+xp+xo)^ 



Hence, the function y{x) in ( 2.26[ ) is given by 

4o-e2(4tp3+(x-t)p+xo), 



y{x) 



'P_ 

-1 _|_ g4(4tp3+(x-t)p+xo) ' 



We c an no w co mpute M according to (2.26) and then find the one-soliton ui{x,t) 
from (2.20). By (2.16), the second column of M equals ^3 for k € DiDD2. Letting 
o" = — 1 we find 



Up {x, t) = — 2p^sech^ (4tp^ + {x — t)p + xq) 
( 

[H^{x,t,k)]2 = 



p^scdn? (^Atp^ +{x—t)p+XQ^ 
2k(k+ip) 



V 



(4.2) 
• (4.3) 



The results for a = 1 are obtained from these expressions by shifting xq — )• xq + «vr/2. 



(4.4) 



This gives the following additional class of solitons, see Figure 4.1 

u^{x, t) = 2p^csch^ [Ap^t + p{-t + x) + xo] . 

The solution Up{x, t) is the well-known sech^-soliton, see Figure 4.2, For KdV this 
soliton has positive height, but for equation ( 1.1 ) it has the opposite sign {u~ {x, t) < 
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15 r 



10 - 



-10 -5 

^2 



10 



Figure 4.1 The csch -soliton Up{x,t) for p = —1/2 and xq = —1 at time t = 0. 
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Figure 4.2 The sech^ -soliton u„ {x,t) for p = —1/2 and xq = —1 at time t = 0. 



for all x,t). The solution u^{x,t) is less standard; it is normally not considered, 
because for each time t it has a singularity at 

-Atp^ + tp - xo 

X = . 

p 

However, for the half- line problem Up{x, t) is a perfectly smooth solution for alH > 
whenever p < —1/2 and < 0, in which case the singularity never enters the positive 
half of the x-axis. 



4.1.1 Asymptotic solitons 

Assume T = oo and let u{x, t) be the solution of an IBV problem on the half- 
line. Assume that u{0,t) decays to zero as t — )• oo. An application of the nonlinear 



steepest descent method of jTj to the RH problem of Theorem 2.8 shows that the 
solution u{x, t) will asymptote for large t into a sum of one-solitons traveling at 
constant speeds; these solitons are generated by the poles of a(A;)r(A;). Away from 
these solitons the solution is of dispersive character (cf. [13 )> i-e. 

u{x,t) = Y,Up„{x,t) + o{l), t^oo, ^ = 0(1), (4.5) 

n 

where we have denoted the poles by {ipn}i, — 1/2 < pi < P2 < • • • < Pa < 0. In 



view of the assumption 2.7, the set of poles of a{k)T{k) coincides with the set of zeros 
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of d{k), unless a zero of d{k) happens to coincide with a zero of B{k); therefore we 
usually say that the solitons are generated by the zeros of d{k) even though strictly 
speaking some of these zeros may cancel with zeros of B{k). 

If u{x, t) satisfies a linearizable boundary condition at x = 0, then a similar 
asymptotic analysis of the RH problem for M in Theorem 3.1 shows that formula 



(4.5) holds also in this case with {ipn}i , —1/2 < < P2 < • • • < PA < 0, denoting 



the poles of a{k)T{k) = 1/A(A;). We draw the following important conclusion: 

For a linearizable IB V problem the asymptotic solitons are generated by 
the zeros of A{k). 

Which solitons can be generated asymptotically for an IBV problem on the half- 
line? We first note that — 1/2 < < p2 < • • • < PA < 0. Hence, whereas equation 



( 1.1 ) admits solitons Up{x, t) for any value of p < 0, only those satisfying — | < p < 
can be generated asymptotically for the half-line problem. Moreover, among the 
solitons with — | < p < 0, the solitons of csch^-type are singular unless p = —1/2. 
We conclude that the solitons that can be generated asymptotically are: (a) the 
sech^-solitons u~{x,t) for — ^ < p < and (b) the csch^-solitons u~{x,t) for p = — ^. 
The interpretation of the requirement p > — 1/2 is clear: the soliton u^{x,t) travels 
at speed 1 — 4p^, so that p < — ^ if and only if the soliton travels to the left; on the 
half-line these left-moving solitons cannot appear asymptotically. 

4.1.2 Spectral functions 



Using the expression (4.3) for /i3, we can derive explicit formulas for the associated 
spectral functions. We will consider the case of u~{x,t); the analogous quantities for 
Up{x,t) can be obtained by shifting xq — >• + i7r/2. First note that 

ipsech^ [4tp^ + ix-t)p + xo] ( 1 1\ , „ 

H3[x,t,k) = — k^O, 



in accordance with (2.6). Evaluating formula iA.M for //3 at x = t = 0, we find 



2k'^ - 2iptanh{xo)k + p^sech^{xo) p'^sech'^{xo) 

= WTWk ' ^^'^ = -^k^T^- ^'-'^ 

After computing s{k) and fisix, t, k) for all x and t, we can compute /i2 by employing 



equation (2.31). Then S{k) is found from (2.33). For a given finite value of T, we 
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find the following expressions: 

A{k-T) = g^2 (fc2^p2) -p )cosh(2pr-8p r)A; 

+ {k^ + p^) cosh (-8T/ + 2rp - 4a;o) + 2{k^ - ipsinh (2pr - 8p^T) k^ 

- ipsinh (-STp^ + 2Tp - 2xo) k^ 

- ip smh.{2xo)k^ + 2p'^k'^ + {k'^ + p^ 



ipsui\i{2xQ)k^ + 2p^A:^ + [k^ + p^) cosh(2xo)A:^ 
+ (^^ + P^) cosh (2 (4T/ - + xo)) /c^ - V sinh (-8Tp^ + 2rp - xq) /c 



(4.7) 



^p-sinh(2xo)A:-e2<4fc^+'=)V ' 



p^sech^ (— 4Tp'^ + — xq) sech^ 



4:^,- i^ft;- ^ p- ) I 
- cosh (2 (4r/ - Tp + xo)) /c^ - + ipsinh (-STp^ + 2rp - 2xo) k 



- ie2<4'='+'=)^p sinh(2xo)A: + e^K^^^+fc) V - / 



In the limit T — )■ oo, for k £ Di U D3 and —1/2 < p < 0, this yields 



A(k;oo) 



2A;2 — 2iptanh(xo)A: +p2gech2(xo) 
2k'^ - 2ikp 



B{k- 00) 



p2sech2(xo) 
' 2/c2 - 2ikp ' 



(4.8) 



As A; — )• 00, we have a{k) — )• 1, 6(/i;) — )■ 0, A{k] T) — )• 1, and B{k; T) — )• in agreement 
with the results of subsection |2.4[ We may also verify directly that 



a{k) 



ip 
'2k 



+ 0(1), b{k) 



ip 
2k 



+ 0(1), k^O, 



in accordance with (2.35), and that 



A{k;T) 



+ 0(1), B{k-T) 



if3 



+ 0(1), 



k^O, 



where 



/3 



1 2pT + sinh {2pT - 8p^T) 
'4^ cosh^ {pT - 4p^T) 



in accordance with (2.36). As another consistency check we note that, if T = 00, it 



is straightforward to verify the global relation (2.37) 



The spectral functions a^{k) and b^{k) on the line can be computed from the 
formula 

s^{k)= lim e*(-'=^+('=+^^')*)^V3(x,0,A:). 



Hence 



b^{k) 
a^{k) 





k+ip 



For the problem on the line the solitons are generated by the zeros of a^{k). As 
expected a^(k) has a zero at k = ip corresponding to the soliton u~{x,t). For the 
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half-line problem with T = oo, the solitons are generated by the poles of the function 
r(A;) defined in (2.42), or, equivalently, by the zeros of d(k). A computation using 
dXel ) and ([47| yields 

_ sech^ {-ATp^ + Tp- xp) 
2k{k + ip) 

Letting T — )■ oo, we obtain 



d{k; co) 



k — ip 



[cosh (2 (4r/ -Tp + xq)) k^ + k^ 
+ ipsinh (-8Tp3 + 2Tp - 2xo) k+p"^]. 

keD2, -l/2<p<0. 



k + ip 

As expected, d{k) has a zero at k = ip corresponding to Up{x, t). 



Finally, the definition (|2.42|) together with the above expressions for the spectral 
ctions, 

r(A;; oo) 



functions, yields (we display the result only in the limit T — )• oo) 

{k + ip)p'^ 



{k — ip) (cosh(2xo) A;^ + k"^ — ip sinh(2xo) A: + p^) ' 



-1/2 <p< 0. 

(4.9) 



4.1.3 Linearizable one-solitons 



The one-solitons u^{x,t) satisfy linearizable boundary conditions of type (b) with x 
given as follows: 

^^|-isech^.o, <r = -l, 
I gcsch xo, cr = 1, 

provided that p is given by p = —1/2. These are stationary (time-independent) 
solitons. 

Let us compute r(A;) as defined in (3.6) when a = —1. We first define the function 
fi{k) by (3.4) with x = —1/2 obtained from (4.10). Using the expressions in (4.6) 
for a{k) and b{k) and recalling that z^i(/c) satisfies (3.1), we compute 

Mk)aiu,{k)) - b{u,ik)) 



h{k)h{vi{k)) - a{vi{k)) 



-4cosh(2xo)/c'' - 4r - 2i sinh(2xo)A; - 1. (4.11) 



Using this in (3.5) and (3.6) we find 



and 



A{k) = -4:k{2k + i) cosh2(xo) 
{2k — z)sech^(xo) 



(4.12) 



{2k + i) (8A:2 + 4ztanh(xo)A; + sech^(xo)) ' 

Note that A(fc) has a zero at k = —i/2 as expected]^ More over, note that the 
expression for T{k) in (4.12) coincides with the expression (4.9) for r(A;;oo) in the 
limit p —7- —1/2. This verifies that the linearizable approach gives the correct result 
for the stationary one-solitons u~{x,t). 



In Theorem 3.1 we considered only the generic case of zeros in D2, whereas this zero lies on the 



boundary of D2- However, the theory also extends to this nongeneric situation, see [20| . 
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Remark 4.3 We have established for the stationary soUtons the equahty 

t{k) = hm lim r(A:;r). (4.13) 

p-s>-l/2T-s>oo 

It is important to take the hmit T — )• cxd on the right-hand side before taking the hmit 



p — )• —1/2. Indeed, letting^ — )• —1/2 in (4.7) before letting T — )• oo, yields expressions 



for A and B different than those obtained by taking the limit p — )• —1/2 of A{k; oo) 
and B{k;oo). This is due to the fact that the stationary soliton with p = —1/2 
does not decay along the t-axis, so that the eigenfunction fii is not well-defined when 
T = oo. If the boundary conditions had been linearizable of type (a), the equality 



(4.13) would have followed from equation (3.8). But the stationary solitons satisfy 
linearizable boundary conditions of type (b). 

We get around this problem by first considering the soliton with — l/2<p<0 
for which /ii is well-defined also for T = oo, and then taking the limit p — )• — 1/2. It 



is noteworthy that, as we saw in (4.12), the linearizable approach immediately gives 
the correct result for T. 

Remark 4.4 It can be shown that the only linearizable IBV problems with initial 
profile given by a one-soliton Up{x,0), x > 0, are the stationary solitons. We prove 
this fact for the case of an initial profile given by u~ ; the case of Up is analogous. Let 
go = n~(0,0) and 52 = ^"^^(0,0) be the boundary values of the initial data. Since 
go 7^ 0, linearizable boundary conditions of type (a) are ruled out. In order for a 
linearizable boundary condition of type (b) to exist we must have 



90 = X, 92 = X + 3x^, X constant. 



Since p <Q and 



2 ,2 -4/(cosh(2x0) - 2) 

go = -2p sech xq, 92 = - 



cosh (xO) 



this requires p = —h, which in turn implies that x is given by (4.10). 



Remark 4.5 These stationary solitons may seem trivial, but actually are important 
for the interpretation of the asymptotics of certain solutions. Indeed, whenever a 
IBV problem with a linearizable boundary condition of type (b) is considered, the 
asymptotics of a solution must necessarily include stationary solitons in order to 



maintain the nonzero boundary value at x = 0. Equation (4.10) implies that for each 
nonzero value of x > ~l/2, there are exactly two stationary solitons with value x 
at x = 0. These two solitons are related by reflection in x = (for x = ~l/2 the 
two solitons coincide). For x > they are of csch^-type and only one of them is 
singularity-free on x > 0. For X < they are of sech^-type. Therefore, given a IBV 
problem with linearizable boundary condition of type (b) for some value of x ^ ~l/2, 
we expect the asymptotics to include one of the stationary solitons with this value of 
xj^ It is not obvious what happens for an IBV problem with x < —1/2; perhaps this 
type of boundary condition will generate an infinite sequence of solitons at x = 0, so 
that a stationary state is never reached. 



^The nonlinear nature of the linearizable boundary condition rules out the possibility of the 
stationary asymptotics being accounted for by a sum of smaller stationary solitons. 
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4.2 Two-solitons 



The two-sohtons of equation ( 1.1 ) and their corresponding eigenfunctions are obtained 
foUowing the steps in subsection 2.3.3 with = 4 and = —ki, k^ = —k2- Let 



C2 



Pi -P2 

9j = -kjX + {kj + Akj)t, 
ki = ipi, k2 = ip2, where the parameters satisf; 



ie^^^P2(pi +P2) 
Pi -P2 



Qj = iOj + Xj, 



1,2, 



pi<0, P2<0, xiGMU M + 



X2 £RU 



+ 



TTl 



Then q, (113)12, (/^3)22 ai"e given by the fohowing expressions, where we use the ab- 
breviations ch and sh for cosh and sinh, respectively: 



/(x,t) 



- pI) (ch(292)pf +pI-pI+ pich(2ei)) 
(pich(Gi)ch(e2) - P2sh(Gi)sh(e2))2 



(4.14) 



ikik + ipi){k + ip2)(pich(ei)ch(G2) - P2Sh(ei)sh(e2))2 

(A:ch(2e2)p? - ip2(P2sh(2Gi) + pish(2e2))pi + k {pj - pi) + A;p2ch(2Gi)) 



3)22 



2 



Ak{k + ipi){k + ip2)(pich(ei)ch(G2) - P2sh(Gi)sh(G2)) 
X [4(pich(Gi)ch(G2) -p2sh(Gi)sh(G2))2A;3 

- 4i [pI - pI) ch(G2)sh(Gi)(pich(Qi)ch(G2) - P2sh(Gi)sh(Q2))/c2 
+ p\k + p^k + (pi -P2){pi + P2)ch.{2Q2)p\k 

- P2 {pI + pI) sh(2Gi)sh(2G2)piA; + pich(2Qi) (2ch(2G2>? +pI- pI) k 
+ mP2 [pI-p\) (p2sh(2Qi) +pish(2G2))]. 

We note that, as /c — )• 0, 

i{Pi- pI) iP2 sinh(2Gi) + pi sinh(2Q2)) / 1 l\ ,^/-|^ 
4yt(picosh(Gi)cosh(G2) -p2sinh(Gi)sinh(G2))2 V-1 "V 



in accordance with (2.6). 



''A real value of the parameter Xj, j — 1, 2, corresponds to a soliton of sech^-type, while if Xj has 
imaginary part ■Ki/2 the corresponding soliton is of csch'^-type. 
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Profile at t = 50 



\d{k; oo)| 





-0.10 -0.05 0.00 0.05 0.10 



Figure 4.6 The two-sohton with xi = X2 = 0, pi = —1/3 and p2 = —1/6 at time 
t = 50 (left) and a contour plot of the function \d{k;oo)\ whose zeros generate the 
asymptotic solitons (right). A dark (light) area indicates a small (large) value of 
\d{k; c<d)\; the x and y axes are labeled by Re A; S (—0.1,0.1) and ImA; G (—0.54,0), 
respectively. The zeros lie at k = —i/3 and k = —i/6 as expected. 



4.2.1 Spectral functions 

For simplicity, we now assume that xi = X2 = 0. Evaluating /X3 at x = i = 0, we find 
the spectral functions 



2{k + ipi){k + ip2y 2{k + ipi){k + ip2)' 

The spectral functions for the problem on the line are given by 

[k + ipi)[k + ip2) 

The zeros of a^{k), as expected, occur at k = ipi and k = ip2. 

The expressions for /i2, S{k), and d(k) are too long to present in general, but we 
can illustrate the results in some special cases. Let pi = —1/3 and p2 = —1/6. Then 



d{k; T) ^ 



T^^2 



2A; (18fc2 _ 9iA: - 1) (3 cosh (^) -Fcosh(f)) 



108 cosh 



+18 cosh ( ^ j k^ + imk^ + 27iisinh f ^ j k^ + 18isinh f ^ j k 



2 



10T\ f 2T\ /10T\ { 2T\ 

+36z sinh ( j A^ + 9z sinh ( — j A^ + 3 cosh ( ) ^ ~ ^^^^ ( ~^ j ^ + 17A 

lorv 



— j +12 {9k^ + A;) cosh f — j +2i sinh f ^ J +^ sinh 
As T — )■ 00, we find, for k £ D2, 

-18A;2 -Qik + l 



27 



d{k; 00) 



-ISk"^ + 9ik + I' 
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with zeros at k = ipi = —i/3 and k = ip2 = —i/6 as expected, see Figure 4.6 Also, 
in the hmit T — )• oo, we have 



^(A;;oo) 



72fc2 + 5 



B{k; oo) 



72/t2 -36iA; + 4' 



72fc2 + 36iA;-4' 
and we can easily verify the global relation 

B{k;oo)a{k) - A{k;oo)b{k) = 0, k £ Di. 

Furthermore, the expressions for the spectral functions have the appropriate behavior 
as A; — )• oo and k ^ 0. 



4.2.2 Linearizable two-soliton 

Retaining the assumption xi = X2 = and letting 

pi = -l/2, P2 = -l/4, 

we obtain the solution 

3 (4 cosh (list- 4x)) + cosh(x) + 3) 



u{x, t) 



4 (cosh {^{t- 4a;)) + 3 cosh (i^(3t + 4x))) 



2 ' 



which satisfies a linearizable boundary condition of type (b) with X = ~|- From the 
expressions for the spectral functions. 



b{k) 
a{k) 



-32fe2+24ifc+4 
32k^+5 
-32fc2+24iA:+4 



we obtain A{k) according to (3.5). The result is 

4 



+ 6ik - 1) . 



As expected, A{k) has two zeros at A; = —i/2 and k = —i/A, see Figure 4.7 



4.3 Rational solitons 



In addition to the solitons with exponential decay, equation (1.1) also admits solitons 



with algebraic decay at infinity. Taking the limit p — )• in the formula (4.4) for the 
csch^-soliton with xq = 0, we find the rational solution 



u[x, t) 



{x-Xo- 



(4.15) 



where we have introduced a new parameter xq in the answer to allow a translational 
shift of the initial data. The corresponding eigenfunction can be found by taking the 
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Profile at t = 50 



01) 




-O.ftL -0 10 -0.05 00 0.05 0.10 



Figure 4.7 The linearizable two-soliton with xi = X2 = 0, pi = —1/2 andp2 = —1/4 
at time t = 50 (left) and a contour plot of the function |A(A;)| whose zeros generate 
the asymptotic solitons (right). The zeros lie at A; = —i/2 and k = — i/4 as expected. 



same limit in (4.3). Similarly, letting xi = X2 



VK 

2 



in the formula (4.14) for the 



two-soliton and taking the limit pi,P2 — ^ 0, we find the rational two-soliton 



u{x, t) 



6{t 



X) 



(t^ -3xf + 3(x^ 



(i3 - 3xt2 3 (a;2 4) t 



Although the rational solitons have singularities when considered on the line, for the 
case of the half-line it is always possible to choose these singularities to lie outside 
the half-line, at least initially. It is therefore natural to consider examples of IBV 
problems with initial data given by the profile of some rational soliton. However, 
these solutions do not have sufficient decay at x = oo for the formalism of section 
[2] to be applicable. For example, the spectral functions a{k) and b{k) for an initial 
profile given by the rational one-soliton (4.15) at time t = with xq < is given by 



b{k)\ 
a{k) 



^ ~^ kxQ 



It is clear that these spectral functions do not have the correct behavior (2.35) as 
/c — )• 0. The study of rational solitons lies beyond the spectral theory defined in 
section [2] and therefore we will not consider these examples any further. 



5 Exponential initial profiles 

Our next class of examples consists of IBV problems whose initial data is given by 



an exponential function. The initial and boundary conditions (1.3)-(1.4) describe an 



IBV problem for equation (1.1). Since the boundary conditions are of the form (b) 



of Theorem |3.1[ they are linearizable. Therefore, if we can find the spectral functions 
a{k) and h{k) from the initial data, then we can obtain the function A(fc) from 



equation (3.5) and the zeros of A(fc) will give us information about the asymptotic 
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behavior of the solution Ur{x,t) for large t. Observe that the initial data Ej.{x) of 
(1.3) grows exponentially as x — >• — oo; hence there is no obvious extension of this 
problem to the line. 



Remark 5.1 The particular form of the initial profiles Er{x) can be motivated as 
follows: The simplest possible exponential profiles with decay as x — t- oo are of 
the form uq{x) = ce"'"^^" for some values of c, a S M, and r > 0. Letting x •= 
uq{0) = ce", the condition that this profile is compatible with a linearizable boundary 
condition at the point x = t = 0, takes the form 

uoxxiO) = X + ^X^ , i-e. ce"r2 = ce" + 3c2e2°. 

Solving this equation for c and discarding the trivial solution c = yields 




i.e. uo{x) = Er{x). 



Hence, the profiles Er{x) defined in ( |1.3[ ) are the simplest exponential initial profiles 
compatible with a linearizable IBV problem. 



5.1 Determination of a{k) and b{k) 

In order to find analytic expressions for the spectral functions a{k) and b{k), we 
will analyze the x-part of the scalar Lax pair for (1.1) presented in Appendix [A| In 
particular, 

^xx = {u + k'^)ip, u{x) = Er{x). (5.1) 



We will then use equation (A.2) to determine the matrix eigenfunction /i3(x,0,A;), 
from which a{k) and b{k) are defined according to 



The general solution of equation ( |5.1[ ) depends on two constants ci , C2 and is given 
by 



99(x, k; ci, C2) =(-l)-^/_2fe ^ 1 (5.2 




where Ia{x) denotes the modified Bessel function of the first kind and T{z) is the 
Euler gamma function. We denote by (p{x, k) the particular eigenfunction with the 
asymptotic behavior ip{x, k) ~ e^^ as x — >• oo. This eigenfunction can be found from 
(5.2) by an appropriate choice of ci and C2. Since the Bessel function /a(x) has the 
asymptotic behavior 



36 



it follows that 99(2;, A;) = ip{x, k; ci, C2), where 



ci 



k/r 



C2 = 0, 



I.e. 



(p{x,k) 




2\/eF^\/r2 



r 1 



2k\ 
r J 



The matrix eigenfunction /i3(x, 0, /c) is obtained from (f{x, k) via equation ( A.2 ). This 
yields 



[//3(x,0,fc)]2 



V 



■--(r2-l)oFi|^;^+2 
6fc(2jA;+r) 

0-^1 



. 2»fc . 

5 r ' 



(r^-l) 



3r2 



where qFi{;-;-) is the confluent hypergeometric function cf. [23j. It can be explic- 
itly verified that the expression for /i3(x,0,A;) has the correct asymptotic behavior, 



fi3{x,0,k) — 7- / as X — 7- 00, and that is satisfies the x-part of the Lax pair (2.1) with 
potential u = Er{x). Evaluating fi^ix, 0, k) at x = yields 



bik) 
a{k) 



1) oFi(;^ + 2;|-3^) 



0-^1 



Qk{ir - 2k) 
2ik_l 1 
r ' 3 3r^ 



(5.3) 



Note that oFi(;ik; z) is analytic as a function of k except for a countable number 
of simple poles at k = in, n = 0,1,2,.... The modified Bessel function luiz) is 
analytic in the whole complex plane as a function of u. In particular, a{k) and b{k) 
are analytic for k & Di U D2 as required. It can be directly verified, in accordance 



with (2.35), that 



a{k) 



+ 0(1), b{k) 



+ 0(1), 



A; ^ 0, 



where 



l/i 



2Vr2-l 



2^/3 



5.2 Asymptotic behavior 



Substituting the expressions for a(A;) and 6(A;) given in (5.3) into (3.5) yields an 
explicit formula for ^(k\ The function A(A;) is defined for k G and its zeros lie 
on the imaginary axis. We denote these zeros by {iVnSx ^ ~l/2 < Pi < P2 < • • • < 
Pat < 0. Numerical computations suggest that the zeros are simple. As described 



in subsection 6.2, the solution Ur{x,t) for large t splits into a sum of one-solitons. 



37 



Ur{x, t) = X]n=i ^Pn(^' + '^(1)- For —1/2 < p„ < 0, the soliton Up^{x, t) must be of 
the form ( [42| ), 

(3^5 = -2p^sech^ (4tp^ + (x - t)p„ + XQn) , 

while if pi = —1/2 it can also be of the form ( |4.4[ ). 

The associated expression for A(k) is rather complicated and not suitable for 
analytical methods. However, numerical methods suggest the following picture, see 
Figure |5.2[ 

• For any value of r, A{k) has a zero at A; = —i/2. 

• For r >> 1, the only zero of A(fc) is k = —i/2. 

• For 0.25 ^ T < 3.33 and r ^ 1, A(A;), in addition to the zero at A; = —i/2, 
has one more zero. This zero enters D2 from below at A; = — z/2 as r decreases 
below 3.33; as r decreases further the zero moves upward along the imaginary 
axis; for r — )• 1 it hits the origin k = and turns around; as r decreases further 
the zero moves downward along the imaginary axis. 

• For 0.15 ^ r < 0.25, A{k), in addition to the zero at A; = —i/2, has two more 
zeros. The new zero enters D2 from above at A; = as r decreases below 0.25; 
as r decreases below 0.15 yet another zero enters D2 from above. 

• As r 4, 0, the number of zeros of A(A;) increases. New zeros keep entering D2 
from above. For r ~ 0.05, in addition to the zero at A: = —i/2, there exist eight 
more zeros. 



Remark 5.3 These observations regarding the zeros of A(A;) can be interpreted in 
terms of asymptotic solitons as follows. The zero at A; = —i/2, which is present for 
all r, generates a stationary soliton. As discussed in subsection 4.1.3, this soliton is 
needed in order to maintain the nonzero boundary value at x = 0. For r > 1, we have 
Xr > 0, so that this soliton is of csch^-type. For < r < 1, we have — 1/3 < Xr < 0, 
so that it is of sech^-type. 

The reason for the appearance of another zero as r decreases below 3.33 is not 
clear, but the following seems a plausible explanation. First assume r > 1. Near the 
boundary, we expect the solution Ur{x,t) to attempt to evolve into the stationary 
soliton of csch^-type with value Er{Q) = Xr at x = 0. This soliton has data p = —1/2 

and xq = — csch~^ ( 1 / 1 V?'^ — 1 ) ; we denote it by u'^{x;r). Comparing the profiles 



Er{x) and u^{x;r), we see that Er{x) lies below u^(x;r), see Figure 5.4 If u{x,t) 
were to evolve from the initial profile Er{x) into the stationary profile u^{x; r), then 
the solution would have to increase near x = 0. In view of the Remark |3.2[ an increase 
in u{x,t) near x = must be accompanied by a decrease of u{x,t) for larger values 
of X. Considering that the shape of the difference Er{x) — u^(x;r) is quite similar 
to that of a sech^-soliton, it is natural to expect that in the process of raising the 
solution from Er{x) to u~^{x;r) near the boundary, another one-soliton u~{x,t) is 
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Figure 5.2 For the values of r specified in tlie first column, the second column 
shows the initial profile Ur{x,0) = Er{x) for the linearizable initial-boundary value 
problem (1.4) (the x and y axes are labeled by x and Er{x), respectively). The 
third column displays contour plots of |A(A;)| (a dark (light) area indicates a small 
(large) value of |A(/c)|; the x and y axes are labeled by Re A; G (—0.1, 0.1) and Imk G 
(—0.54,0), respectively). The fourth column shows the number of solitons exhibited 
asymptotically for large t by the solution Ur{x,t) and the corresponding values of p. 
The solitons correspond to the zeros of A(A;). 
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sent off toward x — t- oo. However, since the only right-traveling solitons are the ones 
with — 1/2 < p < 0, only solitons within this range can be generated. For r >> 1 the 
difference in area between Er{x) and u~^{x;r) is very large, while as r ^ 1 it shrinks 



to zero, see Figure 5.5 It is therefore possible for a right-traveling soliton to appear 
asymptotically whenever r falls below the value for which the size of the difference 
Er{x) —u^{x; r) equals that of the largest right-traveling soliton. As r decreases even 
further, we expect the size of the soliton u~{x,t) to decrease at a rate comparable 
to the rate at which Er{x) — u~^(x;r) decreases. In fact, since we expect the sech^- 
soliton to be accompanied by some dispersive wave trains of predominantly positive 
elevation, it is plausible that the size of the generated soliton is actually larger than 
the value suggested by the area difference 

/•oo 

Ar= [Er{x) - u^{x;r))dx. 



For example, for the box example in the following section we will find that an initial 
profile with u{x,Q)dx = —1, generates an asymptotic soliton with area —1.7404. 
In order to test this interpretation we compute 



u„ (x, t)dx = Ap (5.4) 



and 



Hence, 



f 

Jo 



^ . X , ^^-1 r +/ N , V2r2 + 1 

Er[x)dx = , / u [x;r)dx = ;= — 

3r Jo V3 



V6r2 + 3 + 3 ) _ 1 



r r 



Ar = ^ ^ . (5.5) 

Let Pr denote the value of p for the asymptotic soliton generated by Er{x) for 1 < 
r < 3.33. The first soliton appears when r 3.33 and has area 4-^3.33 = — 2. On the 
other hand, the area difference for r ~ 3.33 is Ar ~ —0.77. In Table[T]we compare the 
area difference Ar with the corresponding size of the generated soliton for some values 
of r. In all cases the area of the generated soliton is about twice that of Ar- This is 
in good agreement with the observations for the box profile already mentioned. 



Table 1: Comparing Ar with the size of the generated soliton. 



r 


Ar 


Pr 




3.33 


-0.77 


-0.5 


-2 


2.8 


-0.77 


-0.28 


-1.12 


2 


-0.23 


-0.11 


-0.44 


1.1 


-0.0041 


-0.002 


-0.0080 
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Figure 5.4 The initial profile Er{x) (solid) and the stationary soliton u^{x;r) 
(dashed) for r = 3.33. 
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Remark 5.6 The interpretation of the observations for < r < 1 is similar, except 
that the stationary sohton now is of sech^-type since Er{0) < 0. Denoting this 
stationary sohton by u~{x;r) we find that the area difference Ar = J^{Er{x) — 
u~ {x; r))dx is given by formula ( 5.5 ) just as before. However, this time when the area 
difference increases beyond that of the largest right-traveling soliton, a collection of 
small solitons is generated, since no single soliton can compensate for this difference. 



6 Box-shaped initial profiles 



In this section we consider the IBV problem for equation ( 1.1 ) with initial and bound- 



ary conditions given by (1.5). The initial profile Pi{x) is a box of height h and length 
L positioned with its left edge at x = xq. The boundary conditions are linearizable 
of type (a) and we can proceed in the same way as in the previous section: (1) We 
find the eigenfunction ij,^{x,0,k) by solving the x-part of the Lax pair with the po- 
tential given by Pi{x). (2) We find the spectral functions a{k) and b{k) by evaluating 
^3(x,0, A;) at a; = 0. (3) We find A(A;) from equation ( |3.2[ ). (4) We investigate the 
asymptotic behavior of ui{x,t) by analyzing numerically the zeros of A{k). 

6.1 Determine a{k) and b{k) 
Consider the matrix equation 

V'x - ikasip = Viip, (6.1) 

where is a 2 x 2 matrix valued function and Vi is defined in terms of the initial data 
uj{x, 0) by equation ( |2.2[ ). If we can find a solution ^^{x, k) of this equation with the 
asymptotic behavior ilj{x,k) ~ f,^kxa■i ^ then the eigenfunction ii^{x,^,k) is related to 

'(/'(x, k) by 

ilj{x,k) = ^i^{x,Q,k)e'^'"'\ (6.2) 



For X > Xq + L, equation ( |6.1| ) is simply 

-01 - ika^ii = 0, (6.3) 

so that 



i}{x,k) = e'^'"'\ x>xo + L. 

Letting 

[ip{x,k)]2 = 



the second column of equation (6.1) for x G {xo,xo + L) reads 

{h - 2k^) V'l + /iV2 - "^ikil^i^ = 0, 
{h - 2k^) ^2 + /iV'i + '2ikip2x = 0. 
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We find 



ipi{x,k) =e k (2ze fc k — 2ik — ie * h + ih + 2e * a + 2a 

2ax\ hc2 

ax / 2ax q q 2Q!a: 2ckx 

ip2[x,k)=e k l-2ie k"^ + 2ik^ + ie fe h-ih + 2e * a + 2a 



4q 



C2^ 

4a 



where a = \/hk'^ — k'^ and ci , C2 are two parameters which need to be determined so 
that ipix, k) is continuous ai x = xq + L. A straightforward computation yields 



ci 



C2 



-1 + e I n,a 



4A;2 (A;2 - h) 



le 



k ( 2 ( -1 + e k ) k^ 



e h + h — 2ill + e k ] a] a 



4A;2 (A;2 - h) 



For rr < xq, the potential uj{x, 0) vanishes and equation (6.1 ) is again given by (6.3). 
We deduce that the second column of ip{x, k) has the form 



{x,k)] 



(m 

\a{k)e 



X < Xq, 



(6.4) 



where a{k) and b{k) are two coefficient functions to be determined. As the notation 
suggests, these coefficient functions coincide with the spectral functions a{k) and b{k). 



Indeed, (6.2) and (6.4) imply that the eigenfunction fj.s{x, 0, k) is given for x < by 



[fi3{x,0,k)y. 



a{k) 



(6.5) 



and evaluation at x = shows that [fi3{0,0, k)]2 = {b{k) , a{k))'^ as required. 

By enforcing continuity of ^(x, k) at x = xq, we deduce the following expressions 
for a{k) and b{k): 



a{k) 



h{i cos{k{L + 2xo)) + sin(A;(L + 2xo))) sinh 



2a(A;) 



-ikL 



2a{k) cosh 



a{k)L 



i (h - 2k^) sinh ( 



2a{k) 



(6.6) 



where a{k) = y/hk'^ — /c*. Note that since the initial profile n/(x,0) vanishes for 
X < 0, the spectral functions a{k) and b{k) coincide with the spectral functions a^{k) 
and h^{k) for the problem on the line; this can also be seen directly from (6.5) together 
with the formula 

s^{k)= lim e-*''^''^V3(x,0,A;). 
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/ = {h,L,xo} 



Initial profile /3/(x) 



\A{k)\ 



of solitons 



h = -l 
L = 1 
xo = 0.1 




-0.10 -0.05 



h = 

L = 




1 

fp^ = -0.449 



0.10 -0.05 



0.05 0.10 



Figure 6.1 The function A(A:) for the IBV problem (1.5) with initial profile f3i{x) 
a box of height h = —1 and length L = 1 has one zero when the box lies some 
distance away from the boundary (xo ^ 1-44), but no zeros when the box lies near 
the boundary (0 < xq ^ 1.44). Hence the boundary has a direct effect on the number 
of asymptotic solitons. 
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In accordance with the general theory of section [2| a{k) — t- 1 and b{k) — )■ as A; — t- oo 
in DiU D2, and 

i^/h sinh ( y/hL] i^/h sinh { \/TiL 

+ 

6.2 Asymptotic behavior 



a{k) = ^ '- + 0{l), b{k) = ^ ^ + 0(1), k^O. 



Using the expressions (6.6) for a{k) and b{k), we obtain an exphcit formula for the 



function A(/c) defined by (3.2). The zeros of A(fc) generate the asymptotic solitons. 
Observe that the function A{k) depends on xq via the spectral function b(k), so that 
the position of the box can affect the position and the number of zeros of A{k). On 
the other hand, the asymptotic solitons for the infinite line problem with the same 
box-shaped initial profile are generated by the zeros of a^{k) = a{k). Since a{k) is 
independent of xq, its zero set is independent of the position of the box. 

In order to see the effect of the boundary explicitly, we consider a box of height 
h = —1 and length L = 1, see Figure |6.1[ In this case we numerically find the 
following: 

• a^{k) has exactly one zero at fc — 0.4351i independently of the value of xq. 

• For < ^ 1.44, A(A;) has no zeros. 

• For 1.44 < Xq, A{k) has one zero. This zero enters D2 from below at A; = —i/2 
as Xq increases above xq ~ 1.44. As xq grows larger, the position of the zero 
moves upward along the imaginary axis. For xq = 2 the zero lies at ~ — 0.449i. 
In the limit xq — )• 00, the zero approaches the zero of a^{k). 

These observations can be understood in terms of asymptotic solitons as follows: 



Equation (5.4) implies that the soliton Up{x, t) has total area u~{x, t)dx = 
Ap. For the problem on the line there is one asymptotic soliton of total area 
approximately —4 x (0.4351) = —1.7404. This area is comparable to the area 
of the box which is —1, and the difference can be explained by the presence 
of dispersive waves of predominantly positive height. The independence of the 
value of xq is a direct consequence of translation invariance. 

For the half-line problem, a small value of xq means that the initial profile is a 
box located very near the boundary at x = 0. The boundary condition forces 
the solution to vanish at x = and this has a suppressing effect on the solution 
near the boundary, hence no solitons have the opportunity to form. 

For 1.44 < Xq, the initial profile is a box located some distance away from the 
boundary. The impact of the vanishing boundary condition at x = is weak 
and the box is able to send away a soliton. As xq grows larger, the impact of 
the boundary decreases. In the limit xq — )• 00, its impact is negligible and the 
result on the line is recovered. 
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I = {h,L,xo} 



Initial profile /3i{x 



tj of solitons 




Figure 6.2 For a wide and shallow initial profile (5i{x) the effect of the boundary on 
the asymptotic solitons is very small. The zeros of the function A(A:) for a box-shaped 
initial profile of height h = —0.15 and length L = 10 are affected only very slightly 
by the position xq of the left edge. The contour plot of |A(/i:)| has been weighted by 
an exponential factor e^^^ for clarity. 
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I = {/i,L,xo} Initial profile Pi{x) 1^(^)1 tt of solitons 




Figure 6.4 For a box-shaped initial profile /5i{x) of positive height (h > 0) there are 
no solitons generated asymptotically. For the half-line problem this follows from the 
absence of zeros for A(A:). 



This interpretation is reinforced if we repeat the same experiment with the wide 



and shallow box of height h = —0.15 and length L = 10, see Figure 6.2 In this case 
we numerically find the following: 

• a^{k) has two zeros at ki « — 0.330i and k2 « — O.llOi indepedently of the 
value of xq. 

• In the limit xq — >• oo, on the half-line, the function A(A;) has the same zeros as 
a^{k). For xq = 2 the zeros lie within 0.001 of their limiting values. 

• Changing xq has essentially no effect on positions of the zeros of A{k). For 
Xq = the zeros lie at ki ~ — 0.335z and k2 ~ — 0.0120i. 

Remark 6.3 The interpretations of these observations are similar to those for the 
box with h = —1 and L = 1 above, except that now the box is very wide and shallow 
and since the effect of the boundary is only felt at the left edge, the impact of the 
boundary condition is substantially reduced. 

As a final conclusion drawn from the numerical observations, we mention that 
whenever h > 0, there appear to be no asymptotic solitons, neither on the half-line 



nor on the line, see Figure 6.4 This observation is in accordance with the fact that 
the right-traveling solitons are waves of depression (see equation ( |4.2| )) and cannot 
be generated by a box of positive height. 



A Scalar Lax pair 



Although it is usually more convenient to use the matrix Lax pair presented in (2.1 ), 
equation ( 1.1 ) can also be formulated as the condition of compatibility of the following 
scalar Lax pair: 

J ^Pxx = {u + k'^)^p, 



(ft = Uxip + {4:k^ -2q- l)ip. 



2 . 
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where A; G C is a spectral parameter and is a constant and t, k) is an eigenfunc- 
tion. The scalar and matrix Lax pairs are related as follows. Let u{x, t) be given 



and let (p{x,t,k) be the solution of (A.l) satisfying (p{x,t,k) 



kx 



~ e""" as X — )• oo. Let 
fi3{x,t,k) be the matrix eigenfunction defined in (2.5) satisfying fis{x,t,k) — t- / as 
X —7- oo. Then 



Hz{x,t,k) 



Lp{x, t, ik) - l(px{x, t, ik) ip{x, t, -ik) - l(px{x, t, -ik) 
Lp{x, t, ik) + Ifxix, t, ik) ip{x, t, -ik) + l(px{x, t, -ik) 



J{-kx+(k+4.k^)t) 

(A.2) 



(^3 
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